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The time-domain response of highly conducting targets following a rapidly terminated electromag- 
netic pulse displays three distinct regimes: early, intermediate and late time. The intermediate and 
late times are characterized by a superposition of exponentially decaying eigenmodes. At early time 
an ever increasing number of rapidly decaying modes contribute, with the result that the scattered 
electric field displays a universal t^ 1 ^ 2 power law which emerges from the diffusive decay of a pattern 
of surface currents induced by the pulse. The power law amplitude reflects the surface geometry of 
the target, a property that may prove useful in buried target classification in geophysical remote 
sensing applications. 

PACS numbers: 03.50.De, 41.20.-q, 41.20.Jb 



Remote detection and classification of buried targets 
is a key goal in a number of important environmental 
geophysical applications, such as toxic waste drum, and 
landmine and unexploded ordnance (UXO) remediation 
PJ. A common tool used for detection of highly con- 
ducting metallic targets is the time-domain electromag- 
netic (TDEM) method, in which an inductive coil is used 
to transmit EM pulses into the ground. Following each 
pulse, the voltage V(t) induced by the scattered field is 
detected by a receiver coil 0. The magnitude and life- 
time of the currents induced in the target, and hence 
V(t), increase with its size and conductivity. Standard 
TDEM sensors are capable of resolving anomalies from 
very small (of order 1 gram) metal targets |3j, and are 
therefore well suited to detection of relatively large buried 
conducting bodies such as UXO. However, since TDEM is 
a very low frequency (typically of order 100 Hz) method, 
its spatial resolution (limited by the target depth and the 
sensor diameter) is also very low 0. Therefore, the raw 
signal amplitude and lifetime provide gross measures of 
the target size and conductivity, but give no direct in- 
formation about its geometry and other physical charac- 
teristics that would enable discrimination between, say, 
UXO and similarly sized clutter. 

Lacking direct target geometry signatures in TDEM 
data (analogous to, e.g., optical and radar images of un- 
occluded targets), one seeks indirect measures via more 
detailed analysis of V(t). This signal contains informa- 
tion about both intrinsic (target size, shape, geometry, 
and other physical characteristics) and extrinsic (rela- 
tive target-sensor orientation, transmitter and receiver 
coil geometries, pulse waveform, etc.) properties, and 
the key to discrimination is the extraction of the former 
from the "background" of the latter. We shall show that 
such an analysis divides naturally into early, intermediate 
and late time domains. Intermediate time is character- 
ized by a finite superposition of exponential decays, the 
slowest of which eventually dominates and defines late 
time. Early time is characterized by an essentially infi- 
nite number of exponential decays which superimpose to 
generate a 1/y/t universal power law divergence in V(t) 
[5j. The importance of this latter interval is greatly en- 



hanced for ferrous targets whose response is so slow that 
it may in fact comprise the full measured range of V(t). 

To focus the discussion consider the following model. 
At low frequencies the dielectric function in the ground 
and in the target is dominated by the its imaginary part, 
e = Airier /u) |6j, where cr(x) is the dc conductivity (in 
Gaussian units), and the Maxwell equations may be re- 
duced to a single equation for the vector potential, 
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with magnetic induction B = V x A and gauge cho- 
sen so that the electric field is E = — (l/c)<9fA. The 
transmitter loop is modelled by the source current den- 
sity js(x, t) = I()(t)Cr(x), where the current Iq consists 
of a periodic sequence of rapidly terminated pulses, and 
Cy(x) defines the transmitter loop. The magnetic field 
is H = B//i, where fj,(x) is the (relative) permeabil- 
ity. The conductivity and permeability are separated into 
background (<7(,(x), ^h(x)) and conducting target (<7 c (x), 
/i c (x)) components, where a Cl \i c vanish outside the tar- 
get volume V c , and it is assumed only that o~b/o~ c <C 1. 

Equation Q is a vector diffusion equation with diffu- 
sion constant D — c 2 / Airier. Typical values are D(, = 
8.0 x 10 10 m 2 /s for a nonmagnetic background with re- 
sistivity of 10 Sim; D c = 2.3 x 10 2 m 2 /s for an aluminum 
target with resistivity 2.8 x 10~ 8 Sim; and D c — 4.0 m 2 /s 
for a steel target with relative permeability 200 and resis- 
tivity 8.9 x 10 -8 Qm. EM signal propagation distance d 
in time t may be estimated via d ~ V Dt. Early-time re- 
sults will require that Tb = R 2 / '-D& <C r c = L 2 /D c , where 
R is the distance between the sensor and the target, and 
L c is the latter's linear size: target -sensor propagation 
time should be instantaneous on the time scale of the 
electrodynamics of the target itself Q. The associated 
condition ii<i c \J Db/D c is easily satisfied for centime- 
ter scale targets at tens of meters depth, and is even less 
stringent for larger targets. The off-ramp time, r r , for the 
transmitted pulse is assumed to satisfy r r <C r c , so that 
pulse termination also occurs essentially instantaneously 
on the scale of the target dynamics (no particular relation 
between r r and Tb is required). 
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In order to further elucidate the various time scales 
in the problem, consider the homogeneous version of l|T[l. 
valid between pulses. The general solution takes the form 
of a superposition of exponentially decaying eigenmodes 



A(x,«) = ^A n aW(x) ( 
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in which to marks the beginning of the free decay window, 
and the mode shapes a^™' and decay rates A„ satisfy the 
eigenvalue equation 
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These modes (which may be orthonormalized by noting 
that y/aaS n ^ are eigenfunctions of a self adjoint operator) 
correspond to special current density patterns j^(x) = 
(A„/c)CT(x)a ( - n - ) (x) with decaying amplitude, but time- 
independent spatial structure Q- The spectrum {A„} is 
bounded below, with fundamental decay rate Ai ~ l/r c 
governed by the target size, but unbounded from above 
0, 0, with more rapidly decaying modes having 
spatial structure on ever smaller scales. Since Ob <C cr c 
currents in the background are negligible compared to 
those in the target, and it follows that A n , as well as 
the internal structure of a'") fx), x £ V c , are essentially 
independent of the background |l2j | . In this sense A„ and 
a^") are intrinsic properties of the target. Explicitly, one 
finds 



An = ^In [ aW*(x)-dl 
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where the transmitter coil here is an idealized ID curve 
Ct- If the coil has N T windings then J (t) = N T i (T), 
where io is the actual current. The voltage measured in 
the receiver loop is then 



V(t) = XX 
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v n = 
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in which Cr is the idealized ID receiver loop and Nr is 
the number of windings. 

The excitation coefficients A n and V n depend on both 
intrinsic (eigenmode) and extrinsic (transmitter/receiver 
loop geometry, position, orientation, etc.) information. 
Given only V(t), absent any information regarding the 
measurement geometry, target classification relies en- 
tirely on the extractable subset of decay rates A„. The 
mathematical problem is equivalent to the famous "Can 
you hear the shape of a drum?" (i.e., to what extent 
is the shape of a struck drumhead determined by its fre- 
quency spectrum?), but is practically much more difficult 



because no analogue of the Fourier transform exists for 
directly estimating the A„. In contrast, if detailed mea- 
surement information is available, direct prediction of the 
amplitudes, and hence of the full signal V{t) is possible. 
A classification scheme may then be developed based on 
a search for the target model that directly minimizes the 
difference between the measured and predicted data , 
thus circumventing the (generally unstable) problem of 
direct estimation of {A„} from noisy data. 

The number of substantially excited modes in de- 
pends on r r (50-100 /is in many commercial systems). 
Roughly, the terminating pulse will excite a subset (de- 
pending on the extrinsic parameters) of those modes with 
A n A r = l/r r . The higher order modes will decay 
very rapidly, but still contribute strongly at early time 
t — to = 0{r r ). Realistically, one can hope to accurately 
compute only the first few hundred modes jToL ITlj . If 
the largest computable decay rate A max is smaller than 
A r , then the interval < t — to < 1/A max will not be 
accurately modelled. For ferrous targets it is often the 
case that 1/Ai exceeds the measurement window and the 
response is entirely early time. The remainder of this 
paper is therefore concerned with the development of a 
complementary theory that deals with this interval. By 
combining this theory with the mode analysis, a compre- 
hensive model of the entire time-domain signal emerges. 

The analysis proceeds in three steps. (1) An "initial 
condition" for the free dynamics, consisting of a pattern 
of currents confined to the surface of the target, is com- 
puted. (2) The time-development of this surface current, 
namely its diffusion into the interior of the target, is com- 
puted. (3) Finally, this solution is used to compute the 
external field generated at the sensor. 

Step 1: The rapid quenching of the transmitter cur- 
rent leads to an outgoing EM pulse that scatters off the 
target in a complicated way, but exits the target region 
by some transient time itr — to + Ttr, with r tr = 0(jb). 
The assumption T& <C t c implies that the internal field 
A(x, to — T r ) just prior to the pulse termination, remains 
essentially fixed during the interval — r r < t — to < T tr , 
responding only in a thin shell near the boundary, dV c . 
More precisely, at high frequencies where the target skin 
depth S c = yj2D c /u is much smaller than the scale of 
tangential variation of A, the internal field near the sur- 
face, with local normal n, takes the form 113 



A(x,w) = A"(r,w)e-l z lV-^/- D % 
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in which z the coordinate along n, and r is orthogonal to 
it, and n • A" = 0. Continuity of n x A implies that A" 
is also the tangential component of the external field. In 
the time domain, © becomes 



A(x,i) = A (x,i) + AA(x,t) 



(J) 



AA(x, t) 



I I -z 2 /iD c (t-t') 

dt'A^v,t')- [l 



t-t' y/AirD c (t-fy 



valid for t — to <C t c , demonstrating the diffusion of the 
signal inwards from the surface. The current density is 
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given by the same expression, but with A" replaced by 
ct c E I = — (a c /c)d t A^ . Integrating over z, at t tI there is 
an effective surface current, 



K(r,i tr ) 



to+r tI 



dt' 

2^ 



Hc(t - t') 



TT^t-A-" (r, t') (8) 



confined to a thin shell with width ^/ D c (r r + r tr ). 

Equation (JSJ) provides a rigorous foundation for K, ex- 
pressing it in terms the external field at the boundary, 
but the latter has no simple form and is generally un- 
known. We now describe an alternate procedure for its 
direct computation via a self-consistency argument. At 
time t tI Q is to be solved with js = 0. Since all back- 
ground transients have died out, the adtA term is of rel- 
ative order D C R 2 / D^L 2 compared to the curl term and 
may be dropped. It follows that H = (u \7 x A = — V$ 
is the gradient of a magnetic potential satisfying 



V ■ (MbV$) = 0. 
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The solution $o to this equation must satisfy appropriate 
boundary conditions on dVt, namely n- (/Xf,Hb — /i c H c ) = 
and n x (H b - H c ) = (4tt/c)K In both cases, 

H c = V x A c is obtained from the initial internal field 
A(x, to — t t ) evaluated at the boundary. The first con- 
dition imposes a unique solution on $ via the Neumann 
boundary condition 

„ fi . v$ = ^n-H c (r,z = (T), (10) 
Mb 

with formal solution 

$o(x)=/ dV ffJV (x,r')^n-H c (r'), (n) 

JdVa Mb 

where gn is the Neumann green function satisfying —V • 
(/ibVgjv) = <5(x — x') with boundary condition n- Vgw = 
0. The second condition determines K: 



K 



4tt 



n x (V$ + He) 
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Step 2: In order to investigate the subsequent evolu- 
tion of the surface current K we take advantage of the 
rapid variation of the fields near the surface with z. Thus, 
the z-derivatives dominate and to leading order in 
the small parameter e = \/D c (t — t t r)/L 2 one need only 
solve the one-dimensional diffusion equation 



D C 8 2 Z -E A 



d t E = 0, z < 0, 



(13) 



with initial condition E(i tr ) = cr c T 1 K^(z). Here E 1 - = 
E — n(n • E) is the tangential part of E, and /i, a, D 
are treated as constants on either side of the boundary. 
Since the external field varies only on the scales L c , R, to 
leading order one has d z 'E ± (z = + ) = 0. Continuity of 
E^ therefore imposes the Neumann boundary condition 
d z E x (z = 0") = 0. The solution to JT3J) is therefore 

2 e -z 2 /4D c (t-t tI ) 

E(x, t) = E (x, t) + z < 0, (14) 

v c y/4nD c (t - t tr ) 



corresponding to a diffusive Ga ussian sprea d with rapid 
z-dependence is on the scale \jD c (t — £ tr ) = tL c . By 
integrating with respect to time, and enforcing the con- 
dition that A should approach the background solution 
Aq for large z/eL c , one obtains 



AA (x,i) = -^K(r) 



AD c (t - 



e -z 2 /4D c (t-t tI ) 

^AnD c {t - t tI ) 



|z|erfc 



y/4D c (t - t t 
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where erfc(z) = e~ s ds is the complementary 

error function. Since K is, in fact, spread over a width 
\/D c Ttr, equations (|14f> and (|15fl are accurate only in the 
range r tr <C t — to <C t c where the precise microscopic 
structure J7J has been washed out by the diffusion kernel. 

Step 3: Equation l|14|) evaluated at z = pro- 
vides the necessary boundary condition for evaluation 
of the external field to leading order in e. Note that 
E(z = 0~) ~ K^/4/x c /<7 c c 2 (t — ttr) diverges, and conti- 
nuity of E^ leads one to expect a corresponding diver- 
gence in the external electric field. We exhibit this for- 
mally through a correction A<!> to the magnetic potential. 
Thus, the normal component of the curl of (|15fl leads to 
the boundary value n- B(z = 0~) = n- (B + AB), where 

n AB = 4y/t - t tI h ■ V x (jJI^~K\ . (16) 

involves only derivatives with respect to the tangential 
coordinate r, and is valid even if fi c , a c vary on the scale 
L c . We therefore obtain $ = $0 + A$, with boundary 
condition — /ifcfi • A$ = n • AB, and hence to 0(e), 

A*(x) = ( d 2 r'g N {^v')— n-AB(r'), (17) 

J Mb 

which is proportional to y/t — t tr - The correction to the 
external vector potential is obtained by solving the aux- 
iliary pair of equations 



V x AA = — m&VA$ 
V-(a 6 AA) = 0, 



(18) 



and is clearly also proportional to \Jt — ttr- The elec- 
tric field correction AE = -(l/c)9 t AA oc (t - t tI )~ 1/2 
therefore has the promised square root early time diver- 
gence. Measurements of magnetic field or voltage (via 
the time derivative of the of the integral of the magnetic 
flux through receiver loop area) follow directly from 11711 . 

We end by illustrating the early time behavior using 
exact analytical results for a homogeneous sphere of ra- 
dius a in a homogeneous background [14j. We consider 
also an initial static transmitted field, so that the initial 
magnetic field is everywhere described by a scalar poten- 
tial. The initial solution is a superposition of spherical 
harmonics, 



^init — 1 lm s [ iJ,m f^/ n \l _i_ Jin („ 



(r/a)', 

btt(r/aY + c^ t (a/r) 1 



r < a 
r > a. 



(19) 
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with cj™ 



"init 



1 - Mc/Mb) 5I+T ( l = 1 corre- 
sponds to the standard case of a uniform illumination 

field, leading to a dipole response) . At t tY the r < a solu- 
tion remains the same, while b l ™ vanishes. The bound- 
ary condition (|10|l then leads to <&o m = c l ™Yi m , with 
c[, m = —l\i c j(l + l)^b. From fT^l surface current is, 



K lrn = JLl 1 



I [l c 



47TO ^ 'I + 1 A*t> 
From and l|17|) one then obtains, 



\/^ + T)X im . (20) 



(21) 



// 6 a \ / + 1 fj, b J V 7r 
and the external fields are given by 



AB 



Ira 



ifj, b 



-0z(t)V x 



/+1 



X, 



(22) 



AA im = w ^^w(;) l+1 X !m , (23) 



which each display \ft cusps, while 



AE = — JUL 

c 



I - 2(i - it,) Vr/ 



i+i 



X 



Irn •> 



(24) 



has the \ j\fi divergence. Here X; m = — + Z)] _1 / 2 x x 
VY; m are the vector spherical harmonics |13f ■ The spatial 
decay rate of the signal increases with I, but all harmonics 
have the same universal power law time-dependence. 

The author is indebted to E. M. Lavely for numerous 
discussions. The support of SERDP, through contract 
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